The Ising spin glass in a transverse field. 
Results of two fermionic models. 
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Abstract 

We analyze the long range Ising spin glass in a transverse field 
r by using Grassmann variables in a field theory where the spin 
operators are represented by bilinear combinations of fermionic 
fields. We compare the results of two fermionic models. In the 
four state (4S)-model the diagonal Sf operator has two vanishing 
eigenvalues, that are suppressed by a restraint in the two states 
(2S)-model. Within a replica symmetric theory and in the static 
approximation we obtain similar results for both models. They 
both exhibit a critical temperature T^iT) that decreases when 
r increases, until it reaches a quantum critical point (QCP) at 
the same value of and they are both unstable under replica 
symmetry breaking in the whole spin glass phase. 
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1 Introduction 



The Ising model in a transverse field is widely studied for being the sim- 
plest system of interacting spins with quantum dynamics. The most striking 
feature is that the competition between thermal and quantum fluctuations re- 
duce the critical temperature up to a point when a quantum phase transition 
occurs at T = 0, at a quantum critical point(QCP)[|I|. We will not discuss 
here the extensive literature on results for several versions of the model, but 
we will concentrate instead in the quantum Ising spin glass in a transverse 
field. This is represented by a Hamiltonian in which only one component of 
the spins, say the z-component, interact among themselves with a random 
interaction while a uniform, constant field F is applied in the transverse x- 
direction. The experimental realizations of this model are the LiHoxYi_xF^ 
compounds 0. 

In the calculation of the quantum mechanical partition function special 
tools are needed to deal with the non-commuting operators forming the 
Hamiltonian. The method more currently used in the study of short-range 
0] and infinite range spin glasses in a transverse field is the Trotter- 
Suzuki formula |0, that maps a system of quantum spins in d-dimensions to 
a classical system of spins in (d + l)-dimensions, and it is suited to perform 
numerical studies. Another way of dealing with the non-commutativity of 
quantum mechanical spin operators is to use Feynman's path integral for- 
mulations 1^, 1^ and to introduce time-ordering by means of an imaginary 
time < r < /3, where (3 is the inverse temperature. The work by Bray and 
Moore P] established the basis for recent developments in the theory of the 
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quantum Heisenberg spin glass [|]. 

A still different functional integral formulation consists in using Grass- 
mann variables to write a field theory with an effective action where the 
spin operators in the Hamiltonian are expressed as bilinear combinations of 



fermions^, |lOl- The advantage of the fermionic formulation is that it has 
a natural application to problems in condensed matter theory, where the 
fermion operators represent electrons that also participate in other physi- 
cal processes, like super conduct ivity|]TT], |l2l and the Kondo effect [|l3l. In 



the present paper we use two fermionic models within a grassmannian field 
theory]^ to analyze the long range Ising spin glass in a transverse field. The 
novelty resides in the method, as all previous results rely on the Trotter- 
Suzuki approximation[|]. A criticism to the fermionic formulation may be 
that the spin eigenstates at each site do not belong to one irreducible repre- 
sentation = ±|, but they are labeled instead by the fermionic occupation 
numbers no- = or 1 , giving two more states with = 0. We call this the 
"four states" (4S) model, and despite the presence of these two unwanted 
states the 4S-Ising spin glass model describes a spin glass transition with the 
same characteristics as the Sherrington - Kirkpatrick (SK) model [14| in a 



replica symmetric theory. A way to get rid of the unwanted states was intro- 



duced before by Wiethege and Sherrington for non-random interactions 
and it consists in fixing the occupation number + rijj by means of an in- 
tegral constraint at every site. We refer to this as the "two states" (2S)-Ising 
model. 

In sect. 2 we analyse the 4S-Ising and 2S-Ising spin glass models in a 
transverse field, within the static approximation in a replica symmetric the- 



ory . The static ansatz neglects time fluctuations and may be considered an 
approximation similar to mean field theory. Numerical Monte Carlo solutions 
of Bray and Moore's equations indicate that the static approximation repro- 



duces the correct results at finite temperatures||16||. When F = the static 
approximation reproduces the exact results obtained by other methods, in 



particular for the 2S-Ising spin glass model we recover SK equations [Q. The 
results in both models are very similar; they both exhibit a critical spin glass 
temperature Tc{T) that decreases when the strength F of the transverse field 
increases, until it reaches a quantum critical point (QCP) at F^, Tc(Fc) = 0. 
The value of Fc is the same for both models and the 4S-Ising and 2S-Ising 
models are identical close to the QCP. We obtained for both models that the 
replica symmetric solution is unstable [0 in the whole spin glass phase, in 
agreement with previous results with the Trotter-Suzuki method]^]. We left 
sect. 3 for discussions. 
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2 The model and results 



The Ising spin glass in a transverse field is represented by the Hamiltonian 



H = -J:J^^s!S]-2^J:sf (1) 

ij i 

where the sum is over the N sites of a lattice and Jij is a random coupling 
among all pairs of spins, with gaussian probability distribution: 



P{Jij) = e--'S^/^'-^YiV/167rJ2 . (2) 
The spin operators are represented by auxiliary fermions fields: 

Si = ^[4'\aii + a\^a^^] (3) 

where the a|^(ajo-) are creation (destruction) operators with fermion anticom- 
mutatiom rules and cr =1 orj, indicates the spin projections. The number 
operators rii^ = a\^ai(j = or 1, then Sf in Eq.(|) has two eigenvalues ±^ 
corresponding to riii = 1 — rij^, and two vanishing eigenvalues when Uii = rii^. 

We shall use the Lagrangian path integral formulation in terms of anti- 
commuting Grassmann fields described in previous publications |P, [l^l; so we 



avoid giving repetitious details. We consider two models: the unrestrained, 



four states model that has been used previously |Tl], |T2|, |T3[, and also the 
two states model of Wiethege and Sherrington where the number operators 
satisfy the restraint rii-^ + riii = 1, what gives S'f = ±|, at every site ||15| 



The partition function in the 4S-model is given by 



(4) 



while in the restrained model it takes the form: 



^25 = 7;[e-^^n<^K-T + ^ii-i)] 

3 

where (3 — is the inverse temperature. 

By using the integral representation for the Kronecker 5-function: 

5{n,^ + n,^ - 1) = ^ fj c^x,e-^["^T+n,;-i] 



we can express Z45 and Z2S in the compact functional integral form 

1 /"^TT 

3 

where: 



(5) 



(6) 



(7) 



(8) 



and fj^j = for the 4S-model while fij = ixj for the 2S-model. Going to 
Fourier representation we introduce the spinors: 



(9) 



and the Pauli matrices: 

1 

1 



-I 

1 



1 
-1 



(10) 
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to write the spin glass part of the action 



where 



Stm = W]l^{u + n)a^±{u) (12) 



with Matsubara's frequencias uj = {2n + 1)tt and Q = 2rmT . In the static 
approximation, we retain just the term f2 = in the sum over the frequency 
Q. 

The transverse part of the action is given by: 



Ar = EV^]H77'MV^,M (13) 
j 

where the inverse propagator is 

= zu; + jjj + Ta^ (14) 
and the total action can be rebuild as 

= Ar + A^s'g (15) 

where A^gQ is the static component of Eq. ([Tl|) . We are now able to follow the 
standard procedures to get the configurational averaged free energy per site 
by using the replica formalism: 

1 , Z(n) - 1 

lim (16) 



where the configurational averaged, replicated, partition function < >c,a 
Z{n) becomes, after averanging over Jij: 
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a/3 



/oo _ 
n 

a/3 

n{n^/^^^.<^e-'^^"A,(g./^)} (17) 



with the rephca index a = 1, 2, .., n, and 



P'jHY^q^pS]^Sf] (18) 

a/3 



We indicate by 5*1 the static component Sj{il = 0) of Eq.(|T2]). We assume 
a rephca symmetric solution of the saddle point equations: 

Qaf^(3 = q qaa = q + X (19) 

where q is the spin glass order parameter and x is related to the static 
susceptibility by x = Tx- 

The sums over a in the spin part of the action produce again quadratic 
terms that can be linearized by introducting new auxiliary fields, with the 
result 

/oo roo 
DzY[ D^Jja{q,X,f^ja,Z,^a) (20) 

-oo Q, J— oo 

where Dy = --^dye~^^^ and 



Ij^ = J D{^i^^)e - (21) 
with 

G]\uj) = r-^{^) + Ka, (22) 
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K = l3J^qz + l3J^i^ (23) 

The gaussian integral over Grassmann variables is straigthforward[^ |TT], 
T3| , giving the result: 

ln(/,g = H{^^ + ^^3f - A.] (24) 

A„ = [f3J^qz + [3J^i^f + {pVf (25) 

The sum over frequencias can be also easily performed |lT| |1^] and we 
obtain 

Ija = 1 + e^^"^" + e''^"2 cosh (26) 
From Eg. ([T7|) , Eg. (pOD and Eg.(^) we obtain at the saddle point: 

j J — oo Q J—oo Ztt 

j dxja [e'^"' + e''^- + 2 cosh ^ A^] } (27) 



For the four states (4S) model there is no restraint and fija = 0, then 
the integrals over Xja egual unity in Eg. (P7|) . For the restrained two states 
(2S) model we have /ijQ, = ixja from Eg.(|l6D, then the integrals over the 
exponential terms identically vanish in Eg. p7|) . We then obtain for the 
model with 2{p + 1) states, p = or 1: 

f3F, = -{f3Jf[f + 2qx]- / Dz\og[2K,{q,x,z)] (28) 

where 

Kpiq, X,z)=p + r Di cosh (29) 

J—oo 



The saddle point equations for the order parameters are: 



/oo I j-CO j-j^ — BY 

Dz— coshVA + — ^sinhVA} -gp (30) 

-oo Kp J-OO A A 2 

/oo 1 roo h 

-oo ^'k'^^ I-oo ^^7K "'""^ ^^^^ 

where h{^,z) is given in Eg. (p3D . We obtain for the de Almeida-Thouless 
eigenvalue [0 and entropy in both models: 

.2 



AT 



/oo 1 poo 

^ Dz—{Kp J cosh x/A + 



sinh VA]-[J" ^^"^ ^'""^ ^^^^ 



Ai ' '"'-oo WA 



Q /-oo 



5 

K 

roo 



The Landau expansion of the free energy in powers of q gives: 

(3Fp = f3F^ + Bpq^ + Cpq' (34) 
where the coefficients are: 

Sp= [Dp-l][Dp-2] 
Cp = ~[Dp - l\{2{Dp - If + ?,{Dp - 2)} (35) 



and 



/oo , 
D^e' cosh ^2x^(3^ J^e + Z^^r^ (37) 
-oo 
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In Eq.(]37|) we need also: 



x(g = o) = xo = ^^V^p-i (38) 

As Cp < 0, the spin glass phase is characterized by -Bp > 0, giving a max- 
imum instead of a minimum |I^ of the free energy. The critical temperature 
is obtained by solving simultaneously: 

n _ P + Ji{l3c) _^ 

' P+Jom 

The numerical results for the critical temperature Tc{T) and the entropy 
5*0 = S{Tc, r) are shown in Fig. 1 for the 4S-model and Fig. 2 for the 2S- 
model. For large values of the transverse field F the 2S-model and 4S-model 
are undistinguishable. The analytic soluction of Eq.(|5PD when = gives 
the critical value Fc = 2\plJ for both models. When F = 0, the equations 
(30)-(33) for the 2S-model (p=0) reproduce the Sherrington-Kirkpatrick [|14| 
results, while for the 4S-model (p=l), we recover our previous results 0. 

Finally, we comment on the de Almeida-Thouless instability. The exact 
soluction for A"^^ in Eq.(|3^) in both limits, F = and F = Fc, shows that 
A^"^ = at the transition point both for the 2S-model and the 4S-model, 
while numerical results confirm that A"^^(Tc) = for both models on the 
critical line Tc(F). This is a correct result and A.T. acknowledges a flaw in a 
previous publication 
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3 Discussion 



We performed a new study of two quantum Ising spin glass models in a trans- 
verse field by means of a path integral formalism where the spin operators are 
represented by bilinear combinations of fermionic fields. All previous results 
in this problem were obtained with the Trotter-Suzuki approximation 0]. In 
the unrestricted four-states (4S)-model the fermionic representation gives for 
the diagonal S'f-operator two eigenvalues = ±^ and two vanishing eigen- 
values, while in the state (2S)-model the vanishing eigenvalues are suppressed 
by means of an integral constraint. The results in both models were obtained 
with the static approximation and the phase diagram coincides with previ- 
ous results with the Trotter- Suzuki method @]. Regarding the de Almeida- 
Thouless instability!!^, we obtained that the replica symmetric solution is 
unstable in the whole spin glass phase. 

In future work we will apply the fermionic representation of the transverse 
Ising spin glass to problems in condensed mather theory and also the replica 
symmetry breaking in the ordered state will be investigated. 
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5 Figure Captions 

Figl. Critical temperature Tc{T) and entropy on the critical line Sc{Tc, 
for the 4S-model. 

Fig2. Same as Figl. for the 2S-model 
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